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Figure 1 0.— Grid and boundary conditions for anisotropic pin array of reference 7, ref (9). 
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SUMMARY 

Flows through idealized pin arrays were investiagted using a quadrilateral grid finite element model and 
the simplified Ergun model to predict leakage flows and pressure drops in brush seals. The models are in good 
agreement in the laminar region with departures in the laminar-turbulent transition region as defined by the 
simplified Ergun model. No local disturbances in the velocity or pressure fields, symptomatic of turbulence were 
found in the numerical results. The simplified model failed to predict the pressure drop of a 32-pin anisotropic array 
unless the gap is taken as the smaller of the anisotropic gaps. Transitional and anisotropic behavior requires further 
investigation. 


NOMENCLATURE 


A 

flow area without bristles (pins) 

Re 

Reynolds number 

Dp 

hydraulic diameter = 1.5 d 

AP 

pressure drop 

d 

bristle (pin) diameter 

<t> 

bristle pack thickness 

di 

shaft diameter 

w 

mass flow rate 

do 

fence diameter 

V 

velocity 

Go 

mass flux without bristles (pins) 

£ 

porosity 

go 

gap between bristles (pins) 

P 

density 

No 

number of bristles (pins) per cm of circumference 

B 

viscosity (dynamic) 

N e 

number of bristles (pins) in a row (circumferential) 

V 

viscosity (kinematic) 

Nx 

number of bristle (pin) rows 





Subscripts 

s solid 
t total 


INTRODUCTION 

Brush seals are effective, compliant, contact seals. Over their lifetime, these seals are subjected to considerable 
wear, bristle displacement, high pressure drops and thermal loads along with unusual operating conditions and 
hystersis. These limitations were recognized by Fergeson (ref. 1) yet he successfully implemented brush seals as 
replacements for some labyrinth seals in gas turbine engines. Other researchers investigated the geometric effects 
e.g., fence height, clearance (Gorelov et al. (ref. 2)), while others concentrated on understanding the nature of brush 
flows (Braun et al. (ref. 3)) and predicting flows as a function of pressure drop (Chupp(4)). Brush seal flows are 
complex and three dimensional with a variety of patterns recognized, figure 1 . The porous fiber bulk flow model 
(Hendricks et al. (ref. 5)) include the effects of bristle motion and provides a value of direct dynamic stiffness. Brush 
seal rotordyanmics coefficients have been assessed by Childs et al. (ref. 6). Preliminary CFD modeling and valida- 
tion have been completed by Kudriavtsev (refs. 7 and 8) and Athavale (ref. 9). Unfortunately, both these efforts have 
been delayed. 

Kudriavtsev and Braun (ref. 7) validated their CFD code using 2-D arrays of pin-cylinders. The agreement be- 
tween experiment and theory is good, figure 2, including simulation of turbine vane cooling. Expanding the CFD 
code to include flow within the cavity and a porous media model of the brush, brush sealing in a gas turbine engine 
was simulated, figure 3, but not validated (ref. 8). 

In Athavale’ s approach(9), the effect of heating and the 3-D bristle geometry was calculated illustrating that 
bristles isolated by rivering flows can be starved for coolant and produce non uniform heating at the interface, figure 
4. Such heating produces nonuniform wear of the bristles and decreases the effectiveness of the brush seal. 

The extension of the Ergun porous flow model (ref. 10) to brush seal flow data provides reasonably good results 
for gases, reference 11. 

Herein we continue to expand porous flow simulation models for brush sealing by comparing preliminary re- 
sults for flows as computed using CFD modeling. 


BULK FLOW MODELING 

The simplified Ergun model of flows in porous media, reference 11, provides useful dimensionless forms and 
insights into design parameters: 


\\i = APp(1.5d) e 3 /(l-e)(t)Go 

(1) 

= 150 /(Re 1 /(l-e)) + 1.75 


>1 = 1.5 God/|i; Go = p 0 V = w/ A 


A = jt(do -df)/4;Dp = 1.5d 

(2) 

= Vopen / Vtotal = 1 - Vs / Vt 

(3) 


= 1 — 7iNo d 2 /(2(l + do/di)(t)cos(0 + cp)) 
and the upper bound on the thickness and number of rows becomes 

(t) = d Nx = d di No/N 0 (4) 
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where No is the number of bristles per unit length as provided by the manufacturer or by micro-examination of the 
brush interface. 

The values of \|/ and Re^ 1 - e ) are calculated from the data set of Carlile et al. (ref. 12), for gases helium, 
air and carbon-dioxide, and overplotted on the results presented by Ergun (ref. 10) as illustrated in figure 5. 

These results are promising and warrent development of a numerical model to extend the range of validity of a 
porous model. 


NUMERICAL MODELING 

In the ideal case of an array of cylinders in cross flow, symmetry allows considerable simplification and only a 
segement of the array need be modeled and calculated. A 6701 quadrilateral finite element grid was constructed with 
grid concentrations within the gaps. Periodic boundary conditions were applied to the inlet and outlet velocity using 
a slug profile to calculate an exit profile that becomes the inlet velocity, with symmetry along the parting planes and 
no slip conditions at the solid interfaces. The grid model and dimensionless boundary conditions for the FEM-flow 
solver, are illustrated in figure 6. The flow field is considered laminar. The dimensionless parameters are defined as 
(refs. 7 and 8): 


Re 2 = 


( <,a y d \ 


Vd 
go ) 


P = 


^p(vd/go) 2 ^ 


( 5 ) 


where V= w/(p A), go = gl = g2 represents the gap at the inlet and exit boundaries, d is the bristle diameter, p the 
density, and p the dynamic viscosity. In this simulation, the pressure drop across Nx bristle rows is assumed to be 
linear and the porosity can be determined by 


£ = l-Vs/Vt = l-7T/(2V3(l + go/d) 2 ) 


( 6 ) 


Where Vs is the total solid, Vt the total volume, and go, the spacing between the bristles and £ the porosity. 
Selected values of (d/go) are 14, 28, and 35 and characterize typical brush configurations. No effect of pressure dif- 
ferential on the bristle, motion of the interface, or elliptical nature of the flow about the bristle or details of the com- 
plex 3-D flow characterizing a brush are condsidered in this model. 


RESULTS 

In the simulation, the Reynolds number (Re 2 ) was varied over a range of 10 to about 400 for each of the three 
(d/go) values (14, 28, 35) and the pressure drop calculated. These computed pressure drops are given as Table 1, 
and plotted on figure 7 along with the simple Ergun model predictions as a function of maximum gap flow velocity. 
The differences between the results as computed using the ALGOR FEM solver and CFD-ACE FD solver were 
within 1 to 2 percent. 

The difference between the FEM predictions and the simple model become pronounced at maximum gap 
velocities above 25 m/s indicating a departure of the modeling methods. While local circulation or turbulence are the 
most likely contributors, considerable attention was given to the computed flow field and not disturbances were 
found over this range in Reynolds numbers, figure 8; the pressure contours are also quite smooth, figure 9. It is true, 
that the FEM model did not admit turbulence and the first-order upwinding scheme may have excessive artificial 
viscosity to prevent circulation; however prior computation experience, references 7 and 8, does not show effects 
of turbulence for the Reynolds number range of Table 1 and little 2-order upwinding effects. For closely packed 
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cylinders (g/d > 20) viscous effects are very strong due to the walls proximity and flow surface to wetted perimeter 
ratio. Recirculations appear only for Re > 1200 and were corroborated by both ALGOL and CFD-ACE using both 
first and second order convective schemes. 

To illustrate the potential departure from laminar flows, the FEM model results are plotted in terms of the Ergun 
model as figure 5. The three parametric lines (d/go = 14, 28, 35) closely parallel the Blake-Koseny relation repre- 
senting laminar flows in porous media up to Reynolds numbers (Re^ of 15 where the effects of turbulence or 
equivalent disturbances cause an increase in the pressure drop for a specific Reynolds number; this corresponds to a 
maximum computed velocity in the gap of about 20 m/s. These differences remain to be studied. 

In further efforts do determine the relations between the simple Ergun model and computational results, in the 
first case, a set of pin arrays as presented in references 7 and 8 were modeled, figure 10. The array has 3 columns 
of 1 1, 10, and 1 1 pins respectively that are spaced (d + go) both axial and transverse giving an anisotropic gap 
(gl = g2) spacing and flow field. The simplified model was applied, but underpredicted the pressure drop by nearly 
a factor of 3, see figure 5, even though the CFD codes of Athavale (ref. 9) and Kudriavtsev and Braun (refs. 7 and 8) 
both predict the experimental data. In the second case, the transverse (d/gl) was set to 35 and the longitudinal (d/g2) 
was set to 51. The results, Table 1 and figure 5, show a nearly 2:1 departure from the simplified Ergun model. The 
reason for failure of the simple model is being investigated. 


Anisotropic Bristle Spacings 

To this point, using the standard definitions for bristle spacings, correlation of flows in anisotropic bristle 
spaced configurations have been unsuccessful. Redefinition of the void fraction in terms of both the gl (transverse 
or normal to the flow) and g2 (longitudnal or in the flow direction) spacings. 


8 — 1 = (71 / 4) /((l + gl / d)(l + g2 / d) cos(cp)) 
cp = arcsin ((l + gl/d) /(2(l + g2 / d))) 


(7) 


did not correlate the results. 

Due to flow symmetry it was reasoned that flow rates would be controlled by g2, if g2 < gl, and vice versa for 
g2 > gl. However, the control is simply assuming a symmetric geometry letting go be the lesser of gl or g2 with 
porosity defined as equation (6). 

For the case of d/gl = 35 and d/g2 =51, the assumed equivalent symmetric geometry is d/go = 51. The resulting 
locus becomes parallel to that noted on figure 5, except shifted through the point (Re y( 1 - e) = 3.7 and \|/ = 42). The 
estimated pressure drop becomes 5.1 psi with the numerical estimate at 6.1 psi. 

Using this same geometric reduction for the oil data of Braun (1993) (ref. 7 and fig. 5), the estimated pressure 
drop becomes 4. 1 psi and the experimental pressure drop is 4 psi. 


Entrance Effects 

Other considerations include the study of flows within in-line tube banks (Athavale, 1995), which showed that 
for laminar flows up to 7 tubes row were necessary for establishing uniform flows and up to 20 for turbulent flows. 
For laminar flows the pressure coefficient diminished monotonically and has a minimum for the turbulent cases. 
Further, some jetting characterized the flow field where the fluid traveled umimpeded up to seven tube rows. A 
similarity between this finding and those for closely spaced orifices in noted where jetting is strongly dependent on 
thermophysical properties, the number of orifices, and their spacing (Hendricks, 1982). 

These results suggest the complexities associated with flows in similar but geometrically simple configurations 
compared to brush seals, are yet to be understood with any great detail. Further without the simplification of symme- 
try, the computations become excessive. The simplified model certainly has many limitations yet it provides the 
designer with first order values for anticipated flows and pressure drops. 
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CONCLUSIONS 


The prediction of pressure drop across an idealized array of pins for a fixed diameter/spacing ratio show good 
agreement with the simplified Ergun model in the laminar regime. 

Departures that occur are indicative of turbulent effects according to the simple model, but careful investigation 
of the computed velocity and pressure fields appear void of perturbations that are indicative of turbulent behavior. 

The simplified model fails to predict, by a factor of 3, the pressure drops associated with both the numerical and 
the experimental results of a 32 pin anisotropic array and by nearly a factor of 2 for an anisotropic array character- 
ized herein. However, if the geometric configuration were assumed symmetric where go is the smaller of gl or g2 
with porosity defined as equation (6), then a better resolution is provided for the experimental data (Braun, 1993) 
and the numerical results herein. 

The departures of the model due to anisotropic behavior, local circulation and transition to turbulence remain to 
be explored. 
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Table 1. Parametric results for FEM-model and simplified Ergun model. Simplified 
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Figure 1 . — Observed flow patterns in brush seals, (a) Rivering. (b) Jetting, (c) Vortical flow, (d) Lateral 
and parallel flow, (e) End-wall flow, (f) Flow at bristle tips, (g) Flow along bristles. 
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Figure 5. — Dimensionless flow loss versus Reynolds number taken from Ergun (1 0) with: (i) superimposed brush seal data for air, carbon dioxide, and 
helium from Carlile et al. (12). (ii) simple Ergun model and numerical model comparisons (iii) anisotropic array results from references 7 and 9; geometry 
is figure 10. Grid background from Bird, R. B., Stewart, W. E.; and Lightfoot, E. N: Transport Phenomena, John Wiley Press, New York, 1960, p. 200. 



(a) 











(C) 





Figure 6. — Grid for FEM-model. (a) Overall, (b) Upper zone, (c) Lower zone. 



Figure 7. — Pressure drop verses maximum gap flow 
velocity at different packing densities. 
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Figure 8. — Velocity field; simulated idealized brush seal with grid shown fig. 6, Re = 100; d/go = 35. 
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Figure 10. — Grid and boundary conditions for anisotropic pin array of reference 7, ref (9). 
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